Data Scenario and Model Hypothesis

Standard fit report for fits of hierSCAL to DERPA data.

Data Scenario: 1956_allFleets

Model Hypothesis: baseModel

Species: Dover, English, Rock

Stocks: HSHG, QCS, WCVI

Final phase convergence diagnostics

Max Gradient: 0.0274234

Objective Function value: -1350.8439587

Time to fit model: 25.4063833

Full complex

Spawning Biomass

Spawning biomass (red line), catch (grey bars), and scaled biomass indices (coloured points), for all stocks and species. Indices are scaled by the fleet catchability, then by the ratio of spawning biomass to vulnerable biomass.

Figure 1: Spawning biomass (red line), catch (grey bars), and scaled biomass indices (coloured points), for all stocks and species. Indices are scaled by the fleet catchability, then by the ratio of spawning biomass to vulnerable biomass.

Fits to indices

Standardised vulnerable biomass (coloured lines) and the scaled and standardised biomass indices they are fitting to (coloured points), for all stocks and species.

Figure 2: Standardised vulnerable biomass (coloured lines) and the scaled and standardised biomass indices they are fitting to (coloured points), for all stocks and species.

Standardised residuals for model fits to biomass indices (coloured points) and a loess smoother with a 20% confidence interval (coloured lines and grey regions), for all fleets, stocks, and species.

Figure 3: Standardised residuals for model fits to biomass indices (coloured points) and a loess smoother with a 20% confidence interval (coloured lines and grey regions), for all fleets, stocks, and species.

Recruitment

Age-1 recruitments for all species and stocks. Equilibrium unfished recruitment $R_0$ is indicated by the horizontal dashed line.

Figure 4: Age-1 recruitments for all species and stocks. Equilibrium unfished recruitment \(R_0\) is indicated by the horizontal dashed line.

Deviations from expected recruitment for all species and stocks.

Figure 5: Deviations from expected recruitment for all species and stocks.

Stock-recruit curves (solid lines) and modeled recruitments (grey points), gridded over species (columns) and stocks (rows).

Figure 6: Stock-recruit curves (solid lines) and modeled recruitments (grey points), gridded over species (columns) and stocks (rows).

Time-varying catchability

Time series of observation model catchability (log-scale, coloured points) for the commercial fleets and Hecate Strait Multispecies Assemblage survey. Lines show the smoothed trend using a Loess smoother

Figure 7: Time series of observation model catchability (log-scale, coloured points) for the commercial fleets and Hecate Strait Multispecies Assemblage survey. Lines show the smoothed trend using a Loess smoother

Fishing mortality

Estimates of fishing mortality from each fleet (coloured points and lines), gridded over species (columns) and stocks (rows). Fishing mortality rates are found using an iterative Newton-Rhapson solver conditioned on the observed catch.

Figure 8: Estimates of fishing mortality from each fleet (coloured points and lines), gridded over species (columns) and stocks (rows). Fishing mortality rates are found using an iterative Newton-Rhapson solver conditioned on the observed catch.

Modeled removals (blue points) using estimated fishing mortality compared to observed catches (open circles), gridded over species (columns) and stocks (rows).

Figure 9: Modeled removals (blue points) using estimated fishing mortality compared to observed catches (open circles), gridded over species (columns) and stocks (rows).

Growth (Probability of Length-at-age)

Probability curves of length-at-age for males (blue) and females (red). Curves show the probability of each length within an age group, and opacity of the lines increases with age.

Figure 10: Probability curves of length-at-age for males (blue) and females (red). Curves show the probability of each length within an age group, and opacity of the lines increases with age.

Selectivity

Selectivity-at-length for each fleet, gridded over species (columns) and stocks (rows).

Figure 11: Selectivity-at-length for each fleet, gridded over species (columns) and stocks (rows).

Selectivity-at-age for each fleet, for females only, gridded over species (columns) and stocks (rows).

Figure 12: Selectivity-at-age for each fleet, for females only, gridded over species (columns) and stocks (rows).

Reference Points

Yield Curves

Equilibrium yield curves as a function of fishing mortality rates, assuming all fishing mortality comes from the modern trawl fleet.

Figure 13: Equilibrium yield curves as a function of fishing mortality rates, assuming all fishing mortality comes from the modern trawl fleet.

Goldilocks Plots

Re-read Hilborn 2018 for the goldilocks plot

Table

Table 1: Biological reference points for each stock and species.
Species Stock B0 R0 M_m M_f h Bmsy Fmsy Umsy MSY SSB_T D_T Dmsy_T
Dover HSHG 17.8 4.22 0.15 0.1 0.76 6.17 0.18 0.17 0.72 7.56 0.42 1.23
Dover QCS 5.97 1.36 0.15 0.1 0.76 2 0.17 0.15 0.24 1.88 0.31 0.94
Dover WCVI 15.19 5.53 0.19 0.12 0.76 5.06 0.21 0.18 0.76 4.43 0.29 0.88
English HSHG 9.64 5.42 0.18 0.13 0.79 3.32 0.36 0.31 0.72 5.73 0.59 1.73
English QCS 0.54 0.35 0.2 0.15 0.79 0.18 0.35 0.3 0.05 0.36 0.67 2
English WCVI 0.85 0.52 0.19 0.14 0.79 0.28 0.35 0.29 0.07 0.35 0.41 1.25
Rock HSHG 15.83 6.28 0.27 0.18 0.73 5.61 0.26 0.22 0.8 8.55 0.54 1.52
Rock QCS 5.84 1.84 0.24 0.16 0.74 2.02 0.23 0.2 0.31 1.78 0.3 0.88
Rock WCVI 1.77 0.42 0.21 0.14 0.72 0.65 0.26 0.23 0.08 0.75 0.42 1.15

Recruitment correlations

Standardised recruitment deviations. Blue indicates a positive deviation, red indicates a negative deviation, and the area of each circle is proportional to the size of the deviation.

Figure 14: Standardised recruitment deviations. Blue indicates a positive deviation, red indicates a negative deviation, and the area of each circle is proportional to the size of the deviation.

Correlation matrix of recruitment deviations. Blue indicates a positive correlation, and red indicates a negative correlation.

Figure 15: Correlation matrix of recruitment deviations. Blue indicates a positive correlation, and red indicates a negative correlation.

Individual species

UNDER CONSTRUCTION

Compositional data

Growth

Optimisation performance

Phase fit table

Table 2: Optimisation performance of hierSCAL for each phase.
phase objFun maxGrad nPar convCode convMsg time mcmcTime
1 991.2709 0.0004778 580 0 relative convergence (4) 0.460400 NA
2 -300.4909 0.0036795 589 0 relative convergence (4) 1.455817 NA
3 -302.8364 0.0114968 591 0 relative convergence (4) 1.283600 NA
4 -514.1385 0.0055692 825 0 relative convergence (4) 1.437517 NA
5 -519.8698 0.0453074 834 0 relative convergence (4) 2.144217 NA
6 -520.2978 0.0157329 843 0 relative convergence (4) 2.751483 NA
7 -596.3466 0.0269245 846 0 relative convergence (4) 2.948783 NA
8 -705.7436 0.0351604 852 0 relative convergence (4) 2.991250 NA
9 -1225.8134 0.0688755 875 0 relative convergence (4) 3.135433 NA
10 -1328.1828 0.0254070 898 0 relative convergence (4) 3.483000 NA
11 -1350.8440 0.0274234 919 0 relative convergence (4) 3.314883 NA
RE NA NA NA NA NA NA NA